An alternative order parameter for the 4-state Potts model 
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We have investigated the dynamic critical behavior of the two-dimensional 4-state Potts model 
using an alternative order parameter first used by Vanderzande [J. Phys. A: Math. Gen. 20, L549 
(1987)] in the study of the Z(5) model. We have estimated the global persistence exponent 9 g by 
following the time evolution of the probability P(t) that the considered order parameter does not 
change its sign up to time t. We have also obtained the critical exponents 8, z, v, and /3 using 
this alternative definition of the order parameter and our results are in complete agreement with 
available values found in literature. 
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I. INTRODUCTION 

In 1989, Janssen et al. Q and Huse Q pointed out, us- 
ing renormalization group techniques and numerical cal- 
culations, respectively, that there is universality and scal- 
ing behavior even at an early stage of the time evolution 
of dynamic systems without conserved order parameter 
(model A in the terminology of Halperin et al. 3]). 

Since then, a great deal of works on phase transi- 
tions and critical phenomena using Monte Carlo sim- 
ulations in the short-time regime have been published 
and their results are in good agreement with theoretical 
predictions and numerical results found in equilibrium 
EHUHBUHElHEi- In addition, the new ap- 
proach has proven to be useful in determining with good 
precision the dynamic exponent z, as well as the new ex- 
ponent which governs the so-called critical initial slip 
|14|. the anomalous behavior of the magnetization when 
the system is quenched to the critical temperature T c . 

The dynamic scaling relation obtained by Janssen et 
al. for the kth moment of the magnetization, extended 
to systems of finite size, is written as 

M {k \t,T,L,m Q ) = b kfi ^M^{b- z t,b 1 / u T,b' 1 L 1 b Xo mo), 

(1) 

where t is the time evolution, b is an arbitrary spatial 
scaling factor, r = (T—T c )/T c is the reduced temperature 
and L is the linear size of the lattice. The exponents 
(3 and v are as usual the equilibrium critical exponents 
associated respectively with the order parameter and the 
correlation length and xq is related to the exponents z, 
0, j3 and v by the equation 

x = 6z + (3 jv. (2) 
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By setting the scaling factor b = t 1 '* and r = in Eq. 
, we obtain the first moment of the magnetization 

M(t)~mot e , (3) 

where tuq represents the initial magnetization of the sys- 
tem. 

Far from equilibrium, another dynamic critical expo- 
nent was proposed by Majumdar et al. fl5| studying 
the behavior of the global persistence probability P(t) 
that the order parameter has not changed its sign up 
to the time t. At criticality P(t) is expected to decay 
algebraically as 

P(t)~tA (4) 

where 9 g is the global persistence exponent. If the time 
evolution would be a Markovian process, then the expo- 
nent g should obey the equation |l5j] 

e q z=-6z + --^-. (5) 

Z V 

However, as shown in several works 0, 0, 0, 0, 0, 
HI El S |H |2j| the exponent g is an indepen- 
dent critical index closely related to the non-markovian 
character of the process. 

In this paper, we estimate the global persistence ex- 
ponent of the 4-state Potts model employing an alterna- 
tive order parameter first used by Vanderzande |27^ in 
the study of the Z(5) model. Our estimate is in complete 
agreement with the results obtained recently for the Ising 
model with three-spin interactions in one direction and 
for the 4-state Potts model j2|J with its traditional order 
parameter (see Eq. ©)• In addition, comparing with es- 
timates for the Ising and 3-state Potts models [TtJ , there 
are consistently increasing values from q = 2 to q = 4. In 
order to check the validity of that alternative order pa- 
rameter, we also estimate the dynamic critical exponents 
9 and z, along with the static critical exponents v and [3. 
Our results [6 = -0.046(9), z = 2.294(3), v = 0.669(6), 
and j3 = 0.0830(6)] are in complete agreement with pre- 
vious results found in the literature. In Section [H] we 
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describe the model and the order parameter. In Section 
lllll we show the short-time scaling relations and present 
our results. In Section llVI we summarize and conclude. 

II. THE MODEL 

The g-state Potts model [2^, is a generalization of 
the Ising model that preserves the next-nearest neighbour 
interaction, works with only two energies (neighboring 
spins are in the same state or not) but permits to put at 
each site any number of states (0 < q < oo). This model 
encloses a quite number of other problems of statistical 
physics. It undergoes a first-order phase transition when 
q > 4 and a continuous phase transiton for q < 4. Thus 
along with the Ising model, the Potts model is an impor- 
tant laboratory to check new theories and algorithms in 
the study of critical phenomena. Although its exact so- 
lution is not known, several results were obtained during 
the last fifty years ^ HJ EJ H . 

The Hamiltonian of the g-state Potts model is given 

by 

0H = -Kj2 S <rw ( 6 ) 

(id) 

where (3 — 1/ksT and fee is the Boltzmann constant, 
represents nearest-neighbor pairs of lattice sites, 
K is the dimensionless ferromagnetic coupling constant 
and <ii is the spin variable which takes the values <Xj = 
0, • • • , q — 1 on the lattice site i. It is well known that 
the critical point of this model is given by [2^| 

K c = log(l + V?)- (7) 

Usually, the order parameter of this model is given by 

Ml(t) = Z^3T) (l>W - !)) (8) 

where L is the linear size of the lattice and d is the di- 
mension of the system. 

In this paper, we use a different definition for the order 
parameter, first proposed by Vanderzande j2]} studying 
the Z(5) model. It can be written as 



M 2 {t) = 




where the average (• • •) is taken over independent initial 
configurations. 

III. RESULTS 

We performed short-time Monte Carlo simulations to 
obtain the critical exponents for the 4-state Potts model. 



Simulations were carried out for square lattices with 
periodic boundary conditions and dimensions L = 120, 
180 and 240. We also used the lattice sizes L = 20, 30, 40, 
50, 60 and 90 just to estimate the exponent z through the 
scaling collapses for different lattice sizes. The estimates 
for each exponent were obtained from five independent 
bins in the critical temperature. For the exponents 9 g , 
z, (3 and v, each bin consisted of 20000 samples, whereas 
for the exponent 9 we have used 100000 samples. When 
estimating the exponent z through the scaling collapses, 
we used 50000 samples. The error bars are fluctuations 
of the averages obtained from those bins. The dynamic 
evolution of the spins is local and updated by the heat- 
bath algorithm. 

In the following sections we show the results for the 
dynamic and static exponents of the 4-state Potts model. 

A. The dynamic critical exponent 9 g 

First of all we are concerned with the global persistence 
probability P(t). It is defined as the probability that the 
global order parameter has not change its sign up to time 
t. For t = 0, the global persistence probability decays 
algebraically as (l5| 

P{t)~t- e \ (10) 

where 9 g is the global persistence exponent. 

In order to estimate the critical exponent 9 g , a sharp 
preparation of the initial states is demanded in order 
to obtain a precise value for the initial magnetization 
mo <C 1. After obtaining the critical exponent 9 g for 
several values of the initial magnetization Too, the final 
value is achieved from the limit too — > 0. 

In this moment, it is worth to explain how obtain a 
small value of too in Eq. ©. First, each site on the 
lattice is occupied by a spin variable which takes the 
values a = 0, 1, 2 or 3 with equal probability. After, the 
magnetization is measured by using Mi (t) and then, the 
variables in the sites are randomly chosen up to obtain 
a null value for the magnetization. The last procedure 
is to change S sites on the lattice in order to obtain the 
desired initial magnetization. It is given simply by 

m = J2 I 11 ) 

and a value of Too is obtained changing S sites occupied 
by a = 2 or 3 and substituting them by a = 0. 

In Fig. ^we show the behavior of the global persistence 
probability for L — 240 and too = 0.000625 in double- 
log scales, together with the behavior of the exponent 9 g 
for too = 0.005, 0.0025, 0.00125 and 0.000625. In order 
to obtain these initial magnetizations for this lattice, 5 
should corresponds to S — 288, 144, 72, and 36. 

The extrapolated value of 9 g , when too — > and L = 
240 is 

9 g = 0.474(7). (12) 
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FIG. 1: The time evolution of the global persistence prob- 
ability P{t) for a lattice size L = 240 and m = 0.000625. 
The error bars calculated over 5 sets of 20000 samples are 
smaller than the symbols. The inset displays the exponent 8 g 
for four different initial magnetizations, as well as its extrap- 
olated values. 



The extrapolated values of 9 g for L = 120 and L = 
180 are shown in Table [I] We emphasize that sometimes 
the value of the initial magnetizations for the lattice size 
L = 180 are slighty different from those of the L = 240. 
This is necessary in order to obtain an integer value of S. 

TABLE I: Global persistence exponent 9 g . 



L 




120 


0.470(5) 


180 


0.473(6) 


240 


0.474(7) 



The estimates obtained with the three lattice sizes 
show that the finite size effects are less than the sta- 
tistical errors. Therefore, we conclude that the values for 
an infinite lattice are within the error bars of our results 
for L — 240. 



B. The dynamic critical exponent 8 

Another critical exponent found only in the nonequi- 
librium state is the exponent 9 that characterizes the 
anomalous behavior of the order parameter in the short- 
time regime. Formerly, apositive value was always asso- 
ciated to this exponent 0, S IH 111 IH HI and the 
phenomenon was known as critical initial slip. However, 
some models can exhibit negative values for the exponent 
9. This is the case, for instance, of the Baxter- Wu model 
[f| and the tricritical Ising model, anticipated by Janssen 
et al. |39| and numerically confirmed by da Silva et al. 
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9 for the 4-state Potts model using the order parameter 
described in Section ILTl (Eg. ©). 

Usually the exponent 9 has been calculated using Eq. 
© or through the autocorrelation 



A{t) ~ t 6 



(13) 



where d is the dimension of the system. In the present 
work however we estimated the exponent 9 using the time 
correlation of the magnetization |37| 



C(t) = (M(0)M(i)>, 



(14) 



which behaves as t when (M(t = 0)) = 0. The average 
is taken over a set of random initial configurations. Ini- 
tially, this approach had shown to be valid only for mod- 
els which exhibit up-down symmetry |3*7| . Nevertheless, 
it has been demonstrated recently that this approach is 
more general and can include models with other symme- 
tries [43. This approach can thus be used for the q =/= 2 
Potts models. 

When compared to the other two techniques (Eq. © 
and Ea. (|13[l ). this method has at least two advantages. 
It does not demand a careful preparation of the initial 
configurations (mo <C 1) neither a delicate limit uiq — > 0, 
as well as the knowledge in advance of the exponent z 
(see Eq. ()13(l ). which is an order of magnitude greater 
than 9. In this case, a small relative error in z induces a 
large error in 9. 

In Fig. |21 we show the time dependence of the time 
correlation C(t) in double- log scales for the system with 
L = 240. The linear fit of this curve leads to the value 
9 = -0.046(9). 



u 



L = 240 



6 = -0.046(9) 



In this paper we reobtain the dynamic critical exponent 
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FIG. 2: The time correlation of the total magnetization for 
samples with {M(t = 0)) =0. The error bars were calculated 
over 5 sets of 100000 samples. 

For the lattice sizes L = 120 and L = 180 we obtained 
respectively 9 = -0.045(8) and 9 = -0.046(8). These 
results are in good agreement with those found for the 
same model using the order parameter of the Eq. © 

mm. 
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C. The dynamic critical exponent z 

The critical exponent z was estimated independently 
by means of two techniques. We began using mixed ini- 
tial conditions, in order to obtain the function Fi{£), 
given by |4l| 



(M(t))l l0=1 



(15) 



where d is the dimension of the system. This approach 
proved to be very efficient in estimating the exp onent z 
for a great number of models [SSIEEilSSEJ- In this 
technique, for different lattice sizes, the double-log curves 
of F2 versus t fall on the same straight line, without any 
rescaling of time, resulting in more precise estimates for 
z. 

The time evolution of F2 is shown on log-scales in Fig. 
□ for L = 240. 




FIG. 3: The time evolution of F2(t). The error bars are 
smaller than the symbols. Each point represents an average 
over 5 sets of 20000 samples. 



The slope of the straight line gives 
d 



yielding 



0.872(1), 



z = 2.294(3). 



(16) 



(17) 



This result is in good agreement with the values 
z = 2.294(6) recently obtained for the Baxter- Wu model 
PI, z = 2.3(1) for the Ising model with multispin interac- 
tions H3, and z = 2.290(3) for the 4-state Potts model 

For L — 120 we have obtained z = 2.296(5) and for 
L = 180 we have obtained z — 2.295(5) indicating that 
the finite size effects are less than the statistical errors. 



The second technique consists of studying the param- 
eter 



R(T, t, L) 



V top 



sign 



1 51 

bottorn 



(S M t),o-S aim ))) (18) 



introduced by de Oliveira |4J|. In this case, the scaling 
relation for T = T c is given by 0] 



R{T = T c , t, L x ) = R{T = T c , b~% 6L1), 



(19) 



with b = L2/L1. This equation shows that the dynamical 
exponent z can be easily estimated by adjusting the time 
rescaling factor b~ z in order to obtain the best scaling 
collapse of the curves for two different lattice sizes. 

Fig. 0]shows the parameter R as a function of the time 
(full lines), as well as the scaling collapse (open circles) 
between different pairs of lattice for samples with ordered 
initial configurations (mo = 1). 




FIG. 4: R(t) vs t with ordered initial configurations mo = 1. 
The full lines show the behavior of R{t) for lattices L — 20, 
30, 40, 50, 60 and 90 (from the bottom to the top) and the 
corresponding time rescaled curves for lattices L = 30, 40, 50, 
60, and 90 (open circles). The exponent z obtained for each 
collapse is shown in Table [2] The error bars, calculated over 
five sets of 50000 samples, are smaller than the symbols. 

The best values of z, obtained through the % 2 test 0] 
for different scaling collapses are shown in Table ITT1 

Our results obtained for the collapse of R(t) are in good 
agreement with our results arising from Fo[t), as well 
as the results for the 4-state Potts model |41|, and the 
Baxter- Wu model which belongs to the same universality 
class p|. 



D. The static critical exponents v and /3 

With the value of the exponent z in hand, we can es- 
timate the static exponent v taking the derivative of the 



5 



TABLE II: Estimates of the dynamical exponent z for the 
best scaling collapse of R(t). 



L 2 .— > 


U 


z 


30 i — > 


20 


2.27(5) 


40 i — y 


30 


2.28(4) 


50 i — > 


40 


2.28(5) 


60 i — > 


50 


2.29(3) 


90 i — ► 


60 


2.28(3) 



logarithm of the order parameter 

M(t, t) = t- f *l l "M(\,t l l v "T) (20) 

with respect to r in the critical point 

d T \nM(t, r)| T=0 = t 1/l/z 9 T diLF(T')|T'=o- (21) 

This equation follows a power law that does not depend 
on L and the function F(t') is a scaling function which 
modifies the power law at r' 7^ but still in the crit- 
ical domain. In numerical simulations we approximate 
the derivative by a finite difference. Our results were ob- 
tained using finite differences of K c ± S with S = 0.001. 
In Fig. El the power law increase of Eq. H21(l is plotted in 
double-log scales for L = 240. 

le-4c 1 1 1 — 1 — 1 1 1 1 1 1 1 1 — r-a 




le-07 - 



FIG. 5: The time evolution of the derivative d T \nM(t, t)| t= o 
on log-log scales in a dynamic process starting from an ordered 
state (mo = 1). The error bars are smaller than the symbols. 
Each point represents an average over 5 sets of 20000 samples. 

From the slope of the curve we estimate the exponent 
1/vz for the three lattice sizes. Using the exponent z 
calculated previously, we obtain v — 0.670(9) for L = 
120, v = 0.668(6) for L = 180, and v = 0.669(6) for 
L = 240. 

Finally, we evaluate the static exponent (3 following the 
decay of the order parameter in initially ordered samples 
(too = 1). At the critical temperature r = 0, the scal- 
ing law of Eq. (|20[) allows one to obtain f3/vz which in 
turn leads to the exponent f3, using the previous result 



obtained for the product vz. In Fig. we show the 
time evolution of the magnetization in double-log scale 
for L = 240. 




= 0.0541(1) 



0.5 — 1 1 1 — • — 1 

10 100 200 

t 

FIG. 6: The time evolution of the magnetization for initially 
ordered samples (mo = 1). The errorbars calculated over 5 
sets of 20000 samples, are smaller than the symbols. 

A linear fit of this straight line gives the value (3/vz = 
0.0541(1) leading to j3 = 0.0830(6). For L = 120 we 
obtained (3 = 0.0834(7) and for L = 180 we obtained 
(3 = 0.0835(4). 

Our results for v and /3 are in goo d ag reement with the 
exact results v — 2/3 and (3 = 1/12 [3lll32l|. 



E. Anomalous dimension xq 

Finally, we calculate the value of the anomalous di- 
mension xq of the magnetization which is introduced to 
describe the dependence of the scaling behavior of the 
initial conditions. It is related to the exponents 9, z, and 
(3/uhy 

x a =0z + f3/is. (22) 

Table ITTT1 shows the values of xo obtained with the ex- 
ponents estimated all along this paper. 

TABLE III: The exponent xq for the 4-state Potts model. 



L 


Xo 


120 


0.021(21) 


180 


0.019(20) 


240 


0.019(23) 



Our results show that the anomalous dimension of the 
4-state Potts model has a null value whose meaning is 
the presence of the marginal operator, i.e., the operator 
which has the scaling dimension equal to dimensional- 
ity of the system and whose effect is not modified under 
renormalization-group operations. An unlike value was 
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found recently by Arashiro et al. 26] for this model and 
for the Ising model with three-spin interactions in one 
direction. 

IV. DISCUSSION AND CONCLUSIONS 

In this paper we revisited the 4-state Potts model in 
order to obtain the global persistence exponent 9 g using 
an order parameter first proposed by Vanderzande |27| 
in the study of the Z(5) model. The results are in good 
agreement with each other. By using this alternative or- 
der parameter, we have also estimated the dynamic crit- 
ical exponents 9 and z, as well as the well-known statical 
exponents v and (3. The exponent 9 was estimated us- 
ing the time correlation of the magnetization, whereas 



to obtain the exponent z we used the function 
which combines simulations performed with different ini- 
tial conditions and scaling collapse for the parameter R 
introduced by de Oliveira. The statical exponents were 
obtained through the scaling relations for the magneti- 
zation and its derivative with respect to the temperature 
at T c . Our results, when compared with available val- 
ues in literature support the reliability of this new order 
parameter. 
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